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Abstract 



We examine the prospects for determining tan /3 from heavy Higgs scalar 
production in the minimal super symmetric standard model at a future e~^e~ 
collider. Our analysis is independent of assumptions of parameter unification, 
and we consider general radiative corrections in the Higgs sector. Bounds are 
presented for = 500 GeV and 1 TeV, several Higgs masses, a variety 
of integrated luminosities and ^-tagging efficiencies, and in scenarios with 
and without super symmetric decays of the Higgs bosons. We find stringent 
constraints for 3 ^ tan/3 < 10, and, for some scenarios, also interesting 
bounds on high tan (3 through tbH^ production. These bounds imply that 
simple Yukawa unifications may be confirmed or excluded. Implications for 
soft scalar mass determination and top squark parameters are also discussed. 
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I. INTRODUCTION 



Supersymmetry (SUSY) is currently a promising framework for understanding the 
physics of electroweak symmetry breaking, and its discovery at future coUider experiments is 
an exciting possibihty. In addition to elucidating weak scale physics, however, the discovery 
of SUSY may also shed light on the mechanism of SUSY breaking, and may even provide 
our first glimpse of physics at the grand unified theory (GUT) and Planck scales. The pro- 
gram of extrapolating weak scale measurements to such high scales will be an extremely 
challenging one, and its success is certainly not guaranteed. What is likely, however, is that 
such a program will require a detailed understanding of the properties of the weak scale 
supersymmetric particles, or, in other words, a precise determination of the various weak 
scale SUSY parameters. 

Of the many SUSY parameters, tan/5, the ratio of Higgs vacuum expectation values, is 
important for a number of reasons. A measurement of tan (3 allows one to determine Yukawa 
couplings, and thereby confirm or exclude the possibilities of mh-mj. unification and SO(IO)- 
like Yukawa unification [jl],|^ . The parameter tan (3 is also required to determine soft SUSY 
breaking scalar masses from the measured physical sfermion masses. Detailed knowledge of 
soft mass parameters may allow us to distinguish various SUSY breaking mechanisms 0. 
Finally, because the parameter tan/3 enters all (neutralino/chargino, sfermion, and scalar 
Higgs) sectors of SUSY theories, a precise measurement of tan (3 from one sector allows one to 
check SUSY relations in other sectors and improves the bounds on many other parameters. 

In this paper we consider the prospects for determining tan j3 from the production of 
Higgs scalars in the general setting of a supersymmetric model with minimal field content. 
(For recent studies in the framework of GUT scenarios, see Ref. 0.) The branching ratios 
of the heavy Higgs bosons are strongly dependent on tan/5 and may be weakly dependent 
on other SUSY parameters. In contrast, almost all other observables that depend on tan/? 
also depend on many additional unknown SUSY parameters, which weakens one's ability 
to determine tan/? precisely and in a model-independent way. In addition, heavy Higgs 
production results in an excess of multi-6 quark events, which, given the excellent 6-tagging 
efficiency and purity now expected to be available at future colliders, allows it to be distin- 
guished from standard model backgrounds. This sector therefore holds the promise of an 
exceptionally clean and powerful determination of tan/5. We will see that Higgs processes 
may provide strong constraints on moderate and high tan/3, regions which are particularly 
interesting for GUTs and soft scalar mass determination. 

In this study, we will consider the experimental setting of a future e~^e~ collider 0-0. 
Such colliders have been shown to be promising for Higgs discovery and study ||7|-|lC|]. In 
particular, we will consider two stages of the proposed Next Linear Collider (NLC): the first 
with ^/s = 500 GeV and design luminosity 50 fb~^/yr, and the second with ^/s = 1 TeV 
and luminosity 100 - 200 fb~Vy^ HSl- W^e will display results for a variety of integrated 
luminosities, detector parameters, and systematic uncertainties. 

Our analysis is intended to estimate the power of particularly promising processes for 
determining tan/5 in a general setting. We do not restrict our attention to specific models 
by assuming minimal supergravity boundary conditions or other parameter unifications. 
Rather, we analyze a number of scenarios by choosing various kinematically accessible heavy 
Higgs masses and consider scenarios with and without supersymmetric Higgs decays. In 
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addition, we will discuss what improvements can be expected if information from other 
experiments and processes are incorporated. It should be stressed, however, that if SUSY 
is discovered, the analysis should be optimized for the particular SUSY parameters realized 
in nature, and the SUSY parameters will best be determined by a global fit to all data. 

This paper is organized as follows. In Sec. we describe the dependence of Higgs 
processes on tan (3. We explain our general treatment of Higgs sector radiative corrections 
and discuss the relevant cross sections and branching ratios. Higgs scalar production may 
be detected in a number of channels; we define the channels we will use and the cuts used to 
isolate these signals in Sec. |TTT[ We then describe our Monte Carlo simulations, and discuss 
the backgrounds and systematic errors entering our analysis. Sec. |IV| contains our results 
for tan j3 bounds in a variety of scenarios without supersymmetric Higgs decays. The effects 
of such decays are discussed in Sec. |V[ Readers who are primarily interested in our results 
are referred to the (tan /3, tan (3') plots in Sees. |^ and 0. Interesting applications of these 
results are contained in Sees. |V| and |V11| . In Sec. |Vlli| , we briefly compare our results to 
those that may be obtained with other processes and present our conclusions. 



II. tan/3 AND THE SUSY HIGGS SECTOR 
A. Definitions and radiative corrections 



We begin by reviewing the necessary details of the scalar Higgs sector |TI|| . We consider 
a supersymmetric model with minimal Higgs content, that is, with two Higgs superfields 

These couple through the superpotential 

W = XeHiLE + XdHiQD - \uH2QU - fxHiH2 , (2) 

where the A are Yukawa couplings and fi is the supersymmetric Higgs mass parameter. The 
ratio of the two Higgs scalar vacuum expectation values is 

tan/?=|^. (3) 

There are four physical Higgs scalars in the theory: the charged scalar H^, the CP-odd A, 
and the two CP-even scalars h and H. The CP-even mass matrix is in general given by 

m\sm'^ f3 + ml cos"^ (3 -(m^ + m|) sin/5cos/? \ /^n 
I — {m\ + m|) sin f3 cos (3 m\ cos^ f3 + m| sin^ (3 j ' ^ ' 

in the basis {Hi,H2), where AAi^ contains all the radiative corrections. The mixing angle 
diagonalizing this matrix, a, enters in Higgs scalar pair production at e~^e~ colliders through 
the vertices ZAH and ZZh, which are proportional to sin(/5 — a), and the vertices ZAh and 
ZZH, which are proportional to cos(/3 — a). In the limit of large rriA, cos(/3 — a) —>■ 0. 
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At tree level, all Higgs scalar masses and interactions are completely determined by tan (3 
and one additional parameter, which is conventionally taken to be the CP-odd mass m^. 
The charged Higgs mass is then given by 

m^± =m\ + , (5) 

and the CP-even masses and a are determined by Eq. with AA4^ = 0. However, the 
relations between Higgs masses and mixings receive radiative corrections. These corrections 



may be large in the CP-even sector |[T^, and are dependent on many additional SUSY 
parameters. Precise measurements of other SUSY parameters, for example, the masses 
and left-right mixing in the top squark sector, may significantly constrain the size of these 
radiative corrections. For most of this study, however, we make the conservative assumption 
that no estimates of their size may be obtained from measurements outside the Higgs sector, 
and we also do not assume that such effects are small. 

Given this framework, the Higgs masses and interactions are all independent quantities, 
and must be determined experimentally. In this study, we will apply the results of previous 
analyses to determine the Higgs masses ||13[, cos^(/5 — a) from a{e~^e~ — > Zh) 0, and 



B{H — > hh) from ffz/P production |]T^. The uncertainties in these measurements will 
be incorporated in our study as systematic errors. Given these measurements, the only 
remaining unknown parameter entering the processes we will study is tan /3. 

By appealing to experimental measurements of Higgs masses and interactions, we do not 
exploit theoretical relations between Higgs masses and mixings to constrain tan j3, and we 
do not restrict the applicability of our analysis to a specific set of parameters or radiative 
corrections. Ultimately, however, we must choose some underlying parameters to study so 
that we may present quantitative results. The choice of parameters is guided by the desire to 
choose scenarios that share qualitative features with a large portion of parameter space. In 
this study, we include the qualitative features of radiative corrections by studying scenarios 
in which such effects are given by the leading mf 1-loop contribution arising from a top- 
stop loop without left-right stop mixing. With this correction, the radiative correction to 



the CP-even Higgs mass squared matrix [12] becomes 



and the Hhh vertex becomes 



^Hhh = 2 sin 2a sin(/5 + a) — cos 2a cos{(3 + a) + 3e— — - cos^ a , (7) 



where 

3 



Stt^ cos^ 9w sin^ P 




(8) 



and mi = ^Jfoijni^ is the geometric mean of the two physical top squark masses. The 
other triple Higgs vertices also receive corrections, but these will not enter our analysis. 
Furthermore, in this case, the tree level relation between mii± and given in Eq. is not 
affected. There are additional possible sources of (typically smaller) radiative corrections. 
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including the bottom squark sector, left-right squark mixing, and the gaugino-Higgsino 



sector |]15[. Our analysis procedure does not assume that these effects are absent, and so 
could be applied to scenarios in which these effects are present as well. Of course, to the 
extent that these effects change the underlying physics, the quantitative results presented 
in the following sections will be modified. 

Finally, aside from radiative corrections to the Higgs masses and field compositions, there 
are corrections to the specific processes we consider. These are 1-loop corrections to the 
cross sections and decay widths, which also may depend on unknown SUSY parameters, 
such as squark masses. We will include these effects as systematic errors, and we discuss 
these errors more fully in Sec. |T|. 

B. Cross sections 

The two body production processes involving Higgs bosons at e^e~ colliders are e^e~ — > 
(7*, Z*) H+H- and e+e" ^ Z* ^ Zh, ZH, Ah, AH. As noted in Sec. |TA|, production of 



ZH and Ah are suppressed by cos^(/? — a). For tjia ~ 200 GeV, this is a large suppression, 
and, although these processes have been included in our simulations, they are statistically 
insignificant for this analysis. 

In this study we will consider two energies for the NLC: y/s = 500 GeV and 1 TeV. 
We choose typical heavy Higgs masses within the kinematically accessible range for each of 
these two energies. For the 500 GeV collider, we consider mH± = 200 GeV, and for the 1 
TeV collider, we consider mu± = 200, 300, and 400 GeV. (Here and in the following, we 
choose to fix mH± rather than the more conventional m^, as the charged Higgs will be seen 
to play the central role in this analysis. The CP-odd masses corresponding to the choices 
above are = 183, 289, and 392 GeV.) 

In Fig. ^ we plot the cross sections for H^H^, AH, and Zh as functions of the center 
of mass energy -y/i for the three values of mH± given above. We have set tan/3 = 5; the 
dependence on tan (3 is very weak for m//± ^ 200 GeV. For fixed tan f3 and , the Higgs 
masses and couplings are determined by including the 1-loop radiative correction given 
in Eqs. and (§) with = 1 TeV. We see that given the NLC design luminosities of 
50 fb^^/yr and 100 - 200 fb~^/yr for the two beam energies, and assuming heavy Higgs 
masses sufficiently below threshold, thousands of events per year will be produced through 
these reactions. 

We will also make use of the three-body processes e+e" — > tbH~ ,ibH^ ||TB|. This process 



takes place through the Feynman diagrams of Fig. g, and is greatly enhanced through 
mfotan/? couplings for large tan/3. The production cross section is plotted in Fig. ^, where 
we have set tan/? = 60. (In calculating the cross sections for this figure, we have required 
Et + Eb > 1.02^5/2 to separate this mode from the two-body production of H^H~ followed 
by H^ tb; the detailed cuts we use in our analysis and experimental simulation will 
be presented below.) We will see that for large tan/3, this mode may be extracted from 
backgrounds. Its sensitivity to large tan/3 will then be useful for placing tan/3 constraints in 
this range. Although we will not study them here, we note that bbA and bbH are similarly 
enhanced for large tan j3, and may also be useful if they can be isolated from backgrounds. 
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C. Branching ratios 



Although the two-body production cross sections are nearly independent of tan/3, the 
heavy Higgs branching ratios are very sensitive to tan (3. The decay width formulas are given 
in the Appendix and the branching ratios for if^, A, and H are plotted as functions of tan (3 
in Figs. ^, ^ and respectively. (Insignificant modes with branching ratios never greater 
than 10~^ are omitted.) In these plots, we fix mii±- The other Higgs masses and mixings 
are then determined as functions of tan/3, including the leading radiative correction of 
Eqs. (H) and (||) with = 1 TeV. Note also that we use rrit = 175 GeV and the running 
mass rrib = 3.2 GeV for the dynamical coupling. We have assumed that all SUSY decay 
modes are suppressed, either kinematically or through mixing angles. SUSY decay modes 
and their effect on our analysis will be discussed in detail in Sec. |V|. 

Several features of the branching ratios are important. Throughout this study, we assume 
mH± > rrit + nib, so the decay — > tb is always open. Given that the decay widths are 
governed by Yukawa couplings, one might expect that this decay mode would be dominant 
for all values of tan (3. In fact, however, if charged Higgses can be pair produced at = 500 
GeV, the phase space suppression for this decay is large, and the branching ratio for H~ — >■ 
rP can be substantial, as may be seen in Fig. ^a. Charged Higgs events with mixed decays, 
H^H~ — >■ tbTUjibfu, will be very useful for determining tan/?. We will see that when the 
branching ratio B{H^ tu) depends strongly on tan/3, roughly for 3 ^ tan/3 ^ 10, we 
will be able to determine tan/? precisely from this channel. For mH± = 400 GeV, as can 
be seen in Fig. ^o, the tb phase space suppression is negligible, and the ratio V{rv) /V{tb) 
approaches ml/ 3ml ~ 0.1 for large tan/?. 

For the neutral Higgs bosons, two features are particularly noteworthy. For mH± = 200 
GeV, ti decays are closed, and we see in Figs. |a and ^ that decays to Higgs and gauge 
bosons are substantial for low tan/?. Such modes will result in 46, 66, and even 86 events 
from AH production, leading to distinct signals. These branching ratios decrease rapidly as 
tan /? increases through moderate values, so again we expect strong determinations of tan (3 
in this region. For mH± = 400 GeV, the tt mode is open and completely dominates for low 
tan/?, as we see in Figs. ^ and ^b. However, as tan/? increases, this mode is suppressed 
by tan^ /3, while the 66 mode is enhanced by tan^ /?. Ignoring phase space suppressions, 
the branching ratios cross roughly at tan/? f» [mt/mhY^'^ ~ 7, so again we expect a strong 
determination of tan (3 when its underlying value is in the middle range. 



III. EXPERIMENTAL SIMULATION 



A. Signals and Cuts 

As seen in Sec. ||, Higgs production leads to events with many 6 quarks. We will exploit 
this feature in conjunction with the excellent 6-tagging efficiency, ~ 60%, that is expected 



to be available at future e+e" colliders |jT^. As many decay channels may be open, Higgs 
production contributes to many types of events. In this section, we first list the eight signal 
channels that we will consider. We then return to them in detail, and explain our motivations 
for choosing them. 

We consider the following channels: 
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1 26 + / + g's + cuts la - le below ("iJ+if"" channel). 

2 26 + / + g's + cuts 2a - 2e below ("t6i7±" channel). 

3 36+1/ (+g's). 

4 36 + 0,2,3,... / (+g's). 

5 46. 

6 46+ 1/ (+g's). 

7 46 + 0, 2, 3, . . . / (+g's) (but not 46). 

8 56 (+ / + g's). 

In this list, "6" and "g" denote hadronic jets with and without a 6 tag, respectively, "/" 
denotes an isolated, energetic e, fj,, or r, and particles enclosed in parentheses are optional. 
In our analysis, we assume that hadronically-decaying r leptons may be identified as leptons, 
ignoring the slight degradation in statistics from multi-prong r decays. Thus, for example, 
channel 6 contains all events with 4 tagged 6's and exactly one e, fi or r, with or without 
untagged jets. We defer the details of the implementation of these requirements in our 
simulations to Sec. [Ill lj| . 



We now discuss each channel in detail. Let us begin by discussing channels 1 and 2, which 
are intended to isolate e+e^ H^H^ thrv and e+e^ — thH"^ tbru, respectively. The 
cross sections for these events are strongly dependent on tan/?: as seen in Fig. ^, the 
branching fractions vary rapidly for moderate tan/3, and we will see below that the cross 
section for tbH^ grows rapidly for large tan /?. These channels therefore allow us to bound 
tan/? if we can reduce backgrounds to low levels. The event shapes of these two channels 
and the largest background e'^e" —>■ ti thlv are given schematically in Fig. |^. We see 
that all three processes have exactly the same final state particle content: 26 + g's, with an 
isolated single lepton. The cross section for ti production is about one order of magnitude 
larger than that of H^H~ production, and hence it is crucial that this background (as well 
as others) be reduced by additional cuts. 

To reduce backgrounds in both channels, we first attempt to reconstruct the W boson 
and top quark. Each candidate event has an untagged hadronic system from decay, 
which we denote "had", and two 6 jets, one of which comes from t decay. Thus, we first 
rescale phad to obtain the correct value of the W boson mass: 

{rwPhadf = m^w ' (9) 

where phad is the four-momentum of the untagged hadronic system measured in the detector, 
and rw is the rescaling factor defined by this equation. In the signal event, rw is close to 1. 
Then, we try to reconstruct the top quark mass from the four-momentum r^^^phad and one of 
the 6 momenta. However, if the 6 from t decay decays semi-leptonically, the neutrino carries 
away some fraction of the momentum. We therefore also define, for each 6 jet, a rescaling 
factor Tf, given by 

{rbPh + rwPhs.Af = ml , (10) 



7 



where pb is the four-momentum of the b jet. For the b coming from the t decay, r?, is almost 
1 in an event without semi-leptonic decay. If the b decays semi-leptonically, the neutrino 
typically carries away about 25% of the total energy in the b rest frame |18[, and may 



become larger than 1. However, even in this case, is less than 1.7 for about 90% of the 
events. We therefore require that, for at least one of the b jets, 0.9 < ri, < 1.7. We then 
identify the b jet most likely to have come from t decay as bi and the other as 62- This is done 
as follows: if only one rb is in the range 0.9 < rb < 1.7, we define the b jet corresponding to 
this Vb as 61, and the other as 62- If both r^'s are in this interval, we identify the b with rb 
closer to 1 as 61, and the other as 62- The untagged hadronic system and bi then form the 
candidate top quark system. 

Given these definitions of rw, fb, 61, and 62, we then impose several kinematic cuts. For 
channel 1, this is quite simple. In the H^H~ pair production event, the untagged hadronic 
system and the 2 6's are all decay products of one if^, and hence their total energy is (in 
principle) equal to the beam energy \\/s. On the other hand, the total hadronic system 
of the ti background has more energy, since the decay products of one top quark alone 
already have energy equal to the beam energy. Thus, we require that the energy of the total 
hadronic system be approximately the beam energy: i-^i— A£^^± < rwEY.aA+'fhiEb^+Eb^ < 
\^/s + AE^±. (The numerical values of the cut parameters A£^^± and those that follow 
depend on ^/s^, they are given in Table |.) Furthermore, a cut on the energy of the candidate 
top quark system also effectively reduces the tt background; we require r^Eh^^ + Vb^Eb-^ < 
^y/s — AE^^. We also impose a cut on the invariant mass of the bottom quark pair m^j^j 
eliminate backgrounds in which the 6 jets arise from Z boson decay {e.g., e~^e~ W^W~ Z): 
^6162 — '^z ~ Am^ or mb^b2 ^ "^z + Am^. Finally, we require that the invariant mass 
of the untagged jets satisfy |mhad — fn-w\ < Ami^/, and that the single lepton have energy 
El > 5 GeV and be isolated with no hadronic activity within a cone of half-angle 20°. In 
summary, for channel 1, we adopt the following cuts: 

la - AE^± < riy^had + n.Eb, + Eb, < i^i + AE+±, with 0.9 < n, < 1.7. 

lb rw/^had + rb,Eb, < - AEf . 

Ic mb^b2 ^ ^z — Aniz or m^jfoj > mz + Am^. 

Id |mhad - "mwl < Amiy. 

le The single lepton must be energetic, Ei > 5 GeV, and isolated, with no hadronic 
activity within a cone of half- angle 20°. 

The situation for channel 2 (the "tfoiJ^" channel) is more complicated, since on-shell 
production of H^H^ pairs is also a background. For this reason, we replace cuts la and 
lb above. In order to make sure that we have only one "on-shell" H^, we impose a cut on 
the energy of the total hadronic system: rviz-Ehad + fb^Eb^ + Eb^ > \\/s + AE^±. This cut 
does not effectively eliminate the ti background, in contrast to cut la, and we must thus 
rely on a cut on the energy of the candidate top quark system to remove the ti background. 
In order to reduce the ti background even if we misidentify the b jet, we modify cut lb and 
instead require that rwE^^^ + rb^Eb^ < \^/s — AEf° and riy-Ehad + '"fc-^fe ^ \\f^~^Ef°. In 
addition, we again use cuts Ic - le. The following is then the complete set of cuts applied 
to channel 2: 
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2a riy^had + n,Eb, + Eb, > l^/s + AE+^, with 0.9 < < 1.7. 

2b ri4.Ehad + n,Eb, <^Vs- A^f and riy^had + n,Eb, < - AEf . 

2c mb^fej < niz — Am^ or mfej^j > + Am^. 

2d |mhad — "^vkI < Amw 

2e The single lepton must be energetic, Ei > 5 GeV, and isolated, with no hadronic 
activity within a cone of half- angle 20°. 

The choice of channels 3 - 8 is motivated by a number of considerations. To exploit the 
6-rich events in Higgs signals, we require many b tags. For channels 3-8, the requirement 
of three to five b tags effectively removes most standard model backgrounds. In fact, AH 
production may result in events with as many as eight b quarks, and so channels requiring 
more than five b tags may also be considered. However, once branching ratios and 6-tagging 
efficiencies are included, such channels suffer from poor statistics and do not improve our 
results. 

In channels 3-7, events with exactly 1 lepton are distinguished from the others. Charged 
Higgs interactions are not plagued by large 1-loop corrections, and therefore provide signals 
that are not subject to large systematic errors. On the other hand, AH events may be 
subject to such uncertainties, in particular in the interaction vertex Hhh. (See Sec. |IIID[ ) 



To avoid contaminating the charged Higgs signal with systematic uncertainties, we would 
like to separate the H^H~ and AH events. This may be achieved for some parameters 
by separating 1/ events, as charged Higgs pair production may produce 11 events through 
H^H^ — > tbtb — > bbbbW~^W~ , where one W decays leptonically and the other hadronically. 
AH events generally do not, unless A,H—>-ti decays are open. 



B. Signal Simulation 

We must now determine the size of the signal in each channel after all branching ratios, 
cuts, and tagging efficiencies are included. The cross sections in channels 3-8 are completely 
determined by branching ratios and the 6-tagging efficiency. In channels 1 and 2, where 
kinematic cuts apply, we must simulate each of the signal events. Events were generated with 
a parton level Monte Carlo event generator, using the helicity amplitude package HELAS 
and phase space sampler BASES [^. For both cases, the spin correlations present in 



the decays were not included. Semi-leptonic b decays were simulated with branching fraction 
24% and energy distribution given in Ref. |jl8|. As stated previously, we use the running 
quark mass values rrit = 175 GeV and rrib = 3.2 GeV in branching fraction calculations. 

Hadronization and detector effects were crudely simulated by smearing the parton 
energies with detector resolutions projected to be available at future e+e" colliders [0: 
^^ad/^ = 40%/\/E © 2%, ai'^-fE = 15%/y/E © 1%, and ap^/p± = 5 x 10"^ © 0.1% for 
muons, with E and p± in GeV. The efficiency of cut 2a is quite sensitive to the hadronic 
calorimeter resolution, but we assume that the characteristics of this calorimeter are well- 
understood. For the purpose of imposing the cuts, an isolated lepton was defined to be one 
with no hadronic parton within a cone of half-angle 20°. In addition, we assume that r 
leptons may be identified as leptons. Finally, initial state radiation was not included. For 
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y/s = 1 TeV, the effects of initial state radiation may be substantial, and the cuts we have 
proposed may require modification. 

The number of signal and background events in each channel is heavily influenced by 
the 6-tagging efficiency and purity. Recently, there have been great improvements in this 
area. In this study, we assume that the probability of tagging a b (c) quark as a 6 quark is 
eb = 60% (ec = 2.6%) [0. (The possibility of tagging light quarks as b quarks is negligible.) 
Uncertainties in these parameters will contribute to our systematic errors (see Sec. |111D[ ). 
In addition, we present results for other 6-tagging efficiencies in Sec. 

With these assumptions, we can now present the expected signal cross sections after cuts 



in each channel. In Figs. ^ - 11, we display these cross sections as functions of a postulated 



tan/3' for fixed underlying values tan/5 = 5 and {y^,mH±) = (500 GeV, 200 GeV), (1 
TeV, 200 GeV), (1 TeV, 300 GeV) and (1 TeV, 400 GeV). These figures are generated as 
follows: we first assume that the underlying values of tan (3 and realized in nature are as 
given above. We then calculate the Higgs masses and compositions, including the radiative 
corrections of Eqs. (||) and (H) with fixed = 1 TeV, as discussed in Sec. [II A| . In particular, 
this determines the values of the physical Higgs masses, cos^(/5 — a), and B{H hh) that 
would presumably be measured. With these quantities then held fixed to their measured 
values, we then consider a hypothetical tan /?', and determine the contributions of the H^H~, 
AH, Zh, and tbH^ signals to the various channels after cuts as a function of tan /?'. The 
contributions to channels 3-8 are determined simply by branching ratios and tagging 
efficiencies. Zh events do not contribute to channels 1 and 2, as they never have exactly one 
lepton. The contributions of the other 3 processes to channels 1 and 2 are determined by 
Monte Carlo simulation. The signal efficiency of the cuts for channel 1 are 48% - 34% (31% 
- 19%, 55% - 62%, and 60% - 64%) for = 500 GeV with m^i = 200 GeV (^/i = 1 TeV 
with mH± = 200, 300, and 400 GeV), as we vary tan/?' from 1 to 100. For channel 2, the 
efficiency is 33% (49%, 79%, and 83%) for = 500 GeV with m^i = 200 GeV {y/s = 1 
TeV with m^i = 200, 300, and 400 GeV) and tan/5' = 60. 

We see in Figs. ^ - |ll] the expected behavior, given the branching ratios shown in Figs. ^, 
H, and |. For mH± = 200 GeV, we see in Figs. H and | that channels 3-8 are large for low 
tan /3, and drop rapidly for increasing values of tan /?. The cross sections in these channels 
are enhanced when the A Zh and H hh branching fractions are large, as is the case 
for low tan (3. In addition, we see that the cross section in channel 1 grows rapidly from 
low tan (3 to high tan /?, as the branching ratio for H~ tv grows. Finally, the cross 
section for tbH^ is virtually non-existent for low tan /?, but grows rapidly for tan /? ^ 20. 
For tan (5 ~ 60, we see that the cross section is large enough to produce tens of events per 
year, allowing a promising determination for high tan /5 if the backgrounds are small. 

Each of the dependencies on tan/3 mentioned above is weakened for larger mH±, as may 
be seen first for mfj± = 300 GeV in Fig. |10[ For mj^± = 400 GeV, the ti decay mode is now 
open. The branching fractions of A Zh and H ^ hh are thus not very large even for low 
tan P, and so the dependence on tan /3 of channels 3-8, though still present, is diluted. The 
cross section for channel 1 also does not rise as much as tan (3 increases, as the cross section 
a{H~^H~ — * tbru) is now no longer enhanced by the phase space suppression of H"^ tb, as 
it was in the = 200 GeV case. Furthermore, because of the smaller branching ratio for 
Tiy, as well as the phase space suppression for the production process e~^e tbH^, 
the cross section for channel 2 at high tan (3 is not as large as in the other cases. 
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C. Backgrounds 



In our analysis, we have included the following standard model backgrounds (cross sec- 
tions in fb before cuts at ^/s = 500 GeV and 1 TeV are given in parentheses): ti (540, 180), 
W+W- (7000, 2700), ZZ (400, 150), tiZ (1.2, 4.7), W+W-Z (40, 60), ZZZ (1.0, 0.85), 
uuti (0.01, 0.55), e+e-tt (0.35, 6.0), uuZZ (0.6, 6.5), e+e-ZZ (1.0, 2.5), e+e-W+W' (250, 
1100), uuW+W- (2.0, 16), euWZ (8.0, 70), and tih (2.0, 3.5). The cross sections for all 
but the last of these processes have been calculated in Ref. p2[, and cross sections for the 



last process may be found in Ref. |]16|. The cross section for tth, in fact, depends on tan/?. 
However, as this cross section is small, the influence of this dependence is rather weak, and 
we treat it as constant. 

The contributions of these backgrounds to channels 3-8 are completely determined by 
their cross sections, branching fractions, and and ec- The estimated background cross 
sections for these channels, as well as those for channels 1 and 2, are given in Table |1[ For 
channels 1 and 2, the kinematic requirements of Sec. [Ill A| must be imposed. By far the 
largest background to these channels before cuts is tt production. To obtain an accurate 
estimate of the contribution of this background after cuts, we have simulated tt events using 
the Monte Carlo program described above, neglecting production-decay spin correlations. 
For \/i = 500 GeV, the signal and background efficiencies as each additional cut is applied 
are given in Tables |T| and |3 for channels 1 and 2. The cuts are seen to be excellent for 
removing the tt background, reducing this background to 0.046 fb (channel 1) and 0.012 
fb (channel 2). For ^/s = 1 TeV, the efficiencies for the tt background are also very small 
(1.7 X 10"'^ for channel 1, and 1.4 x 10~^ for channel 2), and the background is again well 
suppressed (0.045 fb for channel 1, and 0.036 fb for channel 2). 

Given the effectiveness of the cuts for the ti background, we next consider other back- 
grounds. We begin with the backgrounds for y/s = 500 GeV. The processes WWZ luq'qbb 
and euWZ evq'qhh (when the e does not go down the beampipe) are possible backgrounds. 
They are, however, effectively removed by the kinematic cuts, especially Ic and 2c, which are 
designed to eliminate events in which both h quarks originate from the decay of a Z boson. 
After the cuts, these backgrounds are negligible. The background tig thlvg generically 
fails cuts Id and 2d, as the untagged hadronic system consists of the gluon jet and also the 
hadronic decay products of the W , and therefore typically has invariant mass greater than 
niw- For this background to pass the cuts, the gluon jet must mistakenly be included in 
one of the tagged h jets, which greatly reduces the background. Furthermore, even if cuts 
Id and 2d fail to eliminate this background, it may be removed by considering the invariant 
mass of the combined lepton and missing momentum. In the signal, this is ~ , whereas 
in the background it is ~ mw We have not included this background and this cut, but we 
expect the cut to degrade our signal efficiency very slightly. 

Events ti{h, Z) hlvhq'q{h, Z) may contribute to channels 1 and 2 when (/i, Z) decays 
hadronically and only two of the six jets are 6-tagged. However, as in the case of tig, such 
events will fail cuts Id and 2d, as the hadronic system again consists of the hadronic decay 
products of the W and additional hadronic jets. Events tiZ where Z vv may pass the cuts, 
however. In addition, e~^e~ti events may also pass the cuts when both electrons are lost in 
the beampipe. We have not simulated these processes, but rather assume conservatively that 
the kinematic cuts do not further reduce these backgrounds. After including branching ratios 
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and tagging efficiencies, channels 1 and 2 combined receive contributions from tt{Z — > uu) 
and e~^e~tt of 0.038 and 0.056 fb, respectively. We will conservatively assume backgrounds 
of 0.06 fb in each channel. 

For -y/i = 1 TeV, tt{Z vv) may again pass all cuts. After including branching ratios 
and tagging efficiencies, this contributes 0.15 fb. The process e^e~ti also contributes, and 
its cross section is now greatly increased, but now the ti system is often not energetic enough 
to pass cuts la and 2a. The Eti distribution for ^/s = 1.5 TeV is presented in Ref. p2 



From this figure, we may extrapolate to ^/s = 1 TeV to estimate that roughly 40% of the 
e~^e~tt events have > y/s/2. (We have checked that reasonable deviations from this 
value do not significantly change the results presented below.) Including branching ratios 
and tagging efficiencies and the 40% efficiency, this background is 0.38 fb after cuts. Finally, 
the uutt background is now non-negligible, and gives 0.088 fb, conservatively including only 
branching ratios and tagging efficiencies. In summary, for ^/s = 1 TeV, we estimate the 
total backgrounds to channels 1 and 2 combined to be 0.62 fb. We assume backgrounds of 
0.31 fb in each channel. 

AH production may also contribute to channels 1 and 2 when the decay H —>■ W^W~ 
is prominent. Although a signal, this mode has a number of systematic errors that would 
contaminate the H^H~ contribution we are trying to isolate. However, in our Monte Carlo 
simulation, we find that for channel 1, essentially all AH events are eliminated by cut la 
(see Table |T|). For channel 2, only the behavior at large tan/? is critical, and in this region, 
the AH contribution is eliminated by B{H — ^ W^W~) ~ 0. 

Finally, we note that we have not included possible supersymmetric backgrounds. If 
present, such backgrounds certainly require more study. However, it is likely that slepton, 
neutralino, and chargino pair production will be greatly reduced by our demands for multiple 
h tags and the accompanying cuts. Bottom and top squark pair production may be the 
leading SUSY backgrounds, but motivated by the fact that such particles carry SU(3)c 
quantum numbers and so are likely to be heavy, we also do not consider such processes in 
our analysis. 



D. Systematic errors 

In this study, a number of systematic errors must be included. Two important sources 
are uncertainties in the 6-tagging efficiency and the running quark mass m^. In addition, 
however, the determination of tan (3 may be degraded by uncertainties arising from the vir- 
tual effects of other SUSY particles on Higgs processes, which depend on unknown SUSY 
parameters. We will incorporate all such uncertainties in our analysis as systematic uncer- 
tainties, and in this section we describe them and give numerical values for these errors. 
We note, however, that if other measurements are available, these systematic uncertain- 
ties may be greatly reduced. In Sec. |V|, we will consider the beneficial effects that other 
measurements may have on our analysis. 

As 6-tagging plays a central role in our analysis, it is clear that an accurate knowledge 
of the 6-tagging efficiency is important. We have included a systematic uncertainty of ±2% 
for efc [P3|. 



The running h quark mass enters the branching ratio formulas of Eqs. ([Al|) and ( [A2| ) . 
As studied in Ref. [H], a measurement of this parameter that is relatively free of theoretical 
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ambiguities from the bb potential and renormalization group equation evaluations is possible 
using the branching ratio of h. Ref. estimates that B{h ^ tt) may be measured at a 
future e"'"e~ collider to ~ 0.5%, leading to a la error on of 150 MeV, and we therefore 
take TUf, to be in the range ruf, = 3.2 ± 0.15 GeV. 

As noted in Sec. p.1 A| , we do not assume tree level or specific 1-loop relations in the Higgs 
sector, but instead will assume these are all independently measured. The uncertainties in 
these measurements then enter our analysis as systematic errors. We consider the Higgs 
masses and interactions in turn. The h mass will be measured very precisely, and errors 
arising from this measurement are negligible for this study. The charged Higgs mass may 



be measured through its tbtb decay mode. Ref. |T^ finds that given underlying masses of 
mH± = 180 GeV and rrit = 140 GeV, the charged Higgs mass resolution is 16 GeV. This 
measurement is likely to be improved if supplemented by information from the tbru decay 
mode. Here, however, we adopt conservatively the error given in Ref. [0 with appropriate 
rescaling, i.e., we take 6mH± = 16 GeV/^/N, where is the number of tbib charged Higgs 
events multiplied by the efficiency of 3.5% given in Ref. [|1^. Studies have not been conducted 
for the A and H masses. We will assume, however, that their masses may be measured to 
the same accuracy as the charged Higgs. 

The parameter cos^(/5 — a) may be determined by a{e~^e~ — > Zh) to an accuracy of 2% 
0, and we take this as its systematic uncertainty. As given in Eq. (|^), there may also be 



large radiative corrections to the Hhh vertex. In Ref. |]14[, the measurement of this vertex 
through the branching ratio B{H hh) has been considered using the process e^e~ 1/9 H, 
followed by H ^ hh. The cross section a{e~^e~ Huu) is suppressed by cos^(/3 — a), but 
may be significant for low tan /3, the region in which an accurate measurement of B{H — >■ hh) 



is important for this study. For example, in Ref. |14], the cross section for this process 
with tan/3 = 1.5, mn = 200 GeV, and y/s = 500 GeV is shown to be 3 fb, leading to 
hundreds of events per year. Thus, when B{H —>■ hh) is large enough to be important for 
this study, a fairly accurate measurement of its value may be obtained. Without detailed 
studies of backgrounds, it is impossible to determine exactly what bounds may be placed on 
B{H hh); for this study, we simply estimate that B{H — > hh) may be measured with an 
error of 10%, and include this in our systematic errors. It should be noted, however, that if 
the uuH cross section is suppressed, one must turn to AH production, and perform a global 
fit to Xnhh, tan/3 and possibly other parameters. Such a fit is beyond our present analysis. 

Finally, there are radiative corrections to the decay widths and production cross sections. 
Those that depend on standard model parameters are predictable, and so even if large may 
be incorporated in the analysis, once calculated. However, those that depend on unknown 
SUSY parameters are more dangerous. SUSY QCD corrections to the hadronic decay width 
of the charged Higgs have been calculated [^. These studies have shown that the corrections 
may in general be large and of order 40%. However, for 200 GeV ^ mfj± 400 GeV and 
squark masses above 500 GeV, the SUSY QCD correction is reduced to 10 - 20%. SUSY 



QCD corrections to neutral Higgs decay widths have been studied in Ref. p5[, with similar 
results. We therefore include 20% systematic errors for the five decay widths r{H^ —>■ tb), 
T{A bb), T{A tt), T{H bb), and T{H tt). 

The electroweak corrections to the cross section come from diagrams involving squarks, 
neutralinos, charginos, and sleptons. Such corrections have been studied for charged Higgs 
production where effects are found to be typically of order 10%, and may be as large 
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as 25%. However, for a given range of tan (3, the bounds will be determined primarily by 
channels in which the cross sections scale as tan^ f3 or some higher power of tan (3. The 
uncertainty induced in tan (3 is then ^ 5 — 12%, which will be seen to be negligible relative 
to other errors in this study, and we therefore do not include this uncertainty. 



IV. RESULTS 

In this section we present quantitative results for the bounds on tan/3 that may be 
achieved. For now, we assume that SUSY decay modes are absent — such decay modes will 
be considered in Sec. 0. We first discuss how we include statistical and systematic errors 
in our calculations of confidence level contours, and then present results for a variety of 
underlying parameters and experimental assumptions. 

To bound tan (3, we must first select a set of underlying SUSY parameters to determine 
the underlying physics scenario that we hope to constrain. In our framework, as discussed in 
Sec. H, this requires us to choose tan f3 and , and also to fix the radiative corrections. 
This then fixes the number of events that will be observed in each channel. (We assume for 
simplicity that the number of observed events is given by the central value corresponding to 
the underlying parameters.) 

As described above, our analysis is general in that it does not assume any fixed form of the 
Higgs radiative corrections. Thus, to determine tan/? experimentally, we begin by taking 
the Higgs masses, cos^(/5 — a), and B{H — >• hh) to be bounded experimentally through 
the methods described above. Given these measurements, the only remaining unknown 
parameter is tan/?. To determine tan/5, we postulate a hypothetical value tan/?', and 
determine if such a value is consistent with the observed numbers of events in each of our 
eight channel. To quantify this consistency, we define a simple Ax^ variable, 

Ax^-i:i«. (u) 

1=1 ^ Stat + ^ sys 

where i is summed over all channels, and A^, (A^/) is the number of events in channel i 
determined by the underlying (postulated) value tan/? (tan/?'). The quantities al^^^ and 
cTgyg are the statistical and systematic errors for channel z, respectively, and for simplicity, 

we add these in quadrature. The statistical error is cXgtat = The systematic error is 

given by 
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(12) 



where the sum is over the systematic uncertainties in the 12 quantities Pj = e^, rrib, mH±, 
niA, ruH, cos2(/? - a), B{H hh), T{H^ tb), T{A bb), T{A ti), T{H bb), and 
r{H — > tt). The deviations APj are the systematic uncertainties described for each quantity 
Pj in Sec. [HI D| . In the following, we will display Ax^ = 3.84 contours in the (tan /?, tan /?') 



plane, which we will refer to as 95% CL. contours. 

In these plots, the underlying scenario is determined by fixing tan/5, mH±, and = 1 
TeV. Note that the underlying value of rrih therefore varies with tan/?. Of course, once ruh 
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is known, one should consider only scenarios that predict rrih in the experimentally allowed 
range. However, without knowing rrih, we prefer to display results for scenarios with ruh 
given by reasonable radiative corrections (in this case, radiative corrections that may be 
produced by = 1 TeV). 

We first consider the ^/s = 500 GeV collider, and choose a typical kinematically accessible 



charged Higgs mass of = 200 GeV. In Fig. |T^, we display 95% CL. contours in the 
(tan tan /?') plane for four integrated luminosities: 25, 50, 100, and 200 fb~^, or 0.5, 1, 
2, and 4 years at design luminosity. For this plot, we assume = 60%, and have included 
all the systematic errors of Sec. [Ill D| . We expect in this case to bound moderate tan/? 
stringently through the strong dependence of Higgs boson branching fractions on tan j3 in 
this range. In addition, we expect to be able to bound large values of tan/? through the 
process tbH^. These characteristics are evident in Fig. O. As examples, we find that for an 



integrated luminosity of 100 fb ^ and the underlying values of tan /5 listed below, the 95% 
CL. bounds that may be obtained are 



tan/? = 2 
tan/? = 3 
tan/? = 5 
tan/? = 10 
tan (3 = 60 



tan/5' < 2.9 , 
2.5 < tan/3' < 3.6 , 

4.5 < tan/?' < 5.5 , (13) 

7.6 < tan/3' < 30 , 
40 < tan (3' <90 . 



Note that Yukawa coupling constants become non-perturbative below the GUT scale if tan /? 
is too close to 1, or if tan/? is too large (tan/3 ^ 70 — 80). Thus, in much of the parameter 
space that is theoretically interesting, significant constraints on tan (3 may be obtained. 

The above results have interesting implications as tests of Yukawa coupling constant 
unification. For example, mb-rrir unification based on GUTs prefers either large (tan (3 ^ 60) 
or small (tan/3 ^ 2) values of tan/3 |]l|. Furthermore, if we assume simple SO(10)-like 
unification tan/3 is approximately given by rrit/mb ~ 50 — 60, since the Yukawa couplings 
for top and bottom quarks are unified at the GUT scale. As one can see in Fig. |12|, the 
values of tan/3 predicted by these scenarios can be easily distinguished. In addition, the 
stringent constraints on tan (3 available in its moderate range are very useful for soft scalar 
mass determination, as will be seen in Sec. 0. 

If the heavy Higgs bosons are not produced in the first phase of a future e~^e~ collider's 
run, or even if they are, it may be advantageous to increase the beam energy. We consider 
next a ^/s = 1 TeV collider. In Figs. O, IT^, and Il5|, we present results for scenarios with 



this higher beam energy and mH± = 200, 300, and 400 GeV, respectively. We plot contours 
for integrated luminosities of 100, 200, 400, and 800 fb~^. For = 200 GeV, the result 
is dramatically improved over the ^/s = 500 GeV case. This is in many ways a nearly ideal 
scenario for this analysis. As charged Higgs decays to tb are still considerably suppressed 
by phase space, B{H'^ —>■ rv) rises rapidly for increasing tan/3. This, in conjunction with 
the large luminosities that are expected to be available at y/s = 1 TeV, implies that the 
statistical errors in channels 1 and 2 are greatly reduced. We see in this case that the bounds 
are stringent throughout the range of tan /3, and, for example, for an integrated luminosity 
of 100 fb~^ and tan/3 = 60, we may constrain tan/3' to the range 50 < tan/3' < 75. 

For mH± = 300 GeV, the bounds are slightly worse, as B{H'^ tb) is now not highly 
suppressed by phase space, and the number of events in channels 1 and 2 is therefore reduced. 
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In addition, the power of channel 2 is reduced by the great increase in e~^e~tt background 
for a/s = 1 TeV, relative to ^/s = 500 GeV. Improved cuts may be able to reduce this 
background and improve the high tan/3 results. Nevertheless, the bounds are still quite 
strong for moderate tan (3, and interesting determinations of high tan/? are possible for 
large integrated luminosities. Finally, for = 400 GeV, the bounds are again weaker, 
but we are still able to distinguish low, moderate, and high tan /5, and the measurement will 
be useful for soft scalar mass determinations, as we will see below. 

We next consider the dependence of our results on the assumed 6-tagging efficiency. In 
Fig. |TB|, we again consider the case -/s = 500 GeV and mH± = 200 GeV, but plot bounds for 
a fixed luminosity of 100 fb^^ and three 6-tagging efficiencies e?, = 50%, 60%, and 70%. We 
see that the effect of increased e?, is roughly to decrease the integrated luminosity required 
to achieve a certain bound. 

We conclude this section with a discussion of the leading sources of systematic errors in 
the results displayed in Figs. ^ - |1^. As can be seen in the figures, the weakest bounds 
are achieved in the small tan (3 region (tan (3^2 — 3) and the large tan (3 region. In the 
low tan (3 region, and for the = 200 and 300 GeV cases, tan j3 is mainly constrained by 
channel 1 (the ^^H~^H~" channel). The reason for this is that, although channels 3-8 are 
sensitive to variations in tan (3, as may be seen in Figs. ^ - |TI|, these channels require many 
b tags. The uncertainty in thus significantly weakens the bounds from these channels 
and is, in fact, the leading systematic error. However, if the systematic uncertainty in ej, 
is reduced, these channels may improve the constraints. For mH± = 400 GeV, channel 1 
loses its significance since B{H'^ tv) is highly suppressed, and channel 3, which has the 
largest cross section among the inportant processes, becomes the most sensitive one to tan j3. 
Again, the largest systematic error is the uncertainty in e?,. Thus, in the low tan (3 region, 
for all considered, the results can be improved if we can reduce the uncertainty in e^, 
though the statistical error is also non- negligible. 

For large tan/5, tan/5 is constrained only by channel 2 (the "tfoif^" process). In this 
case, the primary source of error is statistical, and for typical luminosities, a reduction of the 
systematic errors does not substantially improve the result. However, for mH± = 200 GeV or 
300 GeV and large tan (3, the systematic uncertainty is not negligible if a high luminosity is 
obtained. In this case, the leading sources of systematic error are mu±, rrify and r{H^ —>■ tb), 
and the results may be noticeably improved if these errors are reduced. In Fig. |1^, we show 
the 95% C.L. contours with all systematic uncertainties omitted for mH± = 300 GeV and 
400 GeV. Comparing these figures with Figs. [T^ and |T3|, we see that the results may be 
improved significantly if the systematic errors are greatly reduced. 

Throughout this analysis, we have assumed that no detailed knowledge of the radiative 
corrections to the Higgs sector may be obtained, and we therefore rely on experimental 
measurements of the various Higgs masses and couplings. However, if these corrections are 
well-understood, for example, through detailed measurements of top squark masses and left- 
right mixing, the results of this analysis may be improved significantly. For example, if the 
radiative corrections are highly constrained, the triple Higgs vertex is essentially a function 
of tan/5 only, and we need not rely on a measurement of B{H hh). There is then a 
strong dependence of the multi-6 cross sections on tan f3. We have analyzed this possibility, 
and find, in particular, that channel 8 is then strongly dependent on tan f3, and this leads to 
marked improvements in the low tan f3 region. This is but one example of how information 
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from other sectors may improve these results. It is clear that other measurements from the 
LHC or NLC may significantly improve the results presented here. 



V. SCENARIOS WITH SUSY DECAY MODES 

Up to this point, we have assumed that Higgs scalars decay only to standard model 
particles. For large Higgs masses, however, decays to supersymmetric particles may be 
allowed p7| , |28[| . In this section, we discuss the effects that decays to sleptons, neutralinos, 



and charginos have on our analysis. Squarks are typically heavy, and so decays to them will 
not be considered here. 

In many models, the right-handed charged sleptons are the lightest sfermions, as their 
masses are not increased by SU(3)c or SU(2)l interactions in the renormalization to low 
energies. We will therefore begin by considering the scenario in which heavy Higgs decays 
to pairs I'^^Ir are open and all other SUSY decays are closed. The scalars A and H"^ may 
decay only to and so their branching fractions are unchanged in the absence of Ir-Il 
mixing. On the other hand, the H boson may decay into 1*rIr pairs through a D-term 
interaction. The important point is that the HI*^Ir vertex is completely fixed by the U(1)y 
gauge coupling constant and tan/3, i.e., it does not depend on additional unknown SUSY 
parameters. In addition, the slepton masses will be very accurately measured at future e+e~ 
colliders |2l|j2^. Therefore, in the case where H —>■ 1*j^Ir is the only relevant SUSY decay 



mode, no new systematic uncertainties enter our analysis. The primary effect of this decay, 
then, is only to decrease the number of the signal events from AH production. 

In Fig. the branching ratio of H —>■ 1*rIr is given by the solid curve for fixed m//± = 
300 GeV and three degenerate generations of right-handed sleptons with masses m^-^ = 
100 GeV. We see that the branching fraction never exceeds 0.3 and decreases rapidly for 
increasing tan/3 as the width to h quarks becomes dominant. In particular, for tan/5 ^ 3, 
the range in which our analysis may give stringent bounds, the branching ratio is less than 
0.1. Furthermore, if the decay mode H ^ ti is open, the branching ratio for this SUSY 
decay mode is suppressed even for the low tan /? region. 

If only decays to left-handed slepton pairs are possible, again only the H branching 
fractions are affected. In Fig. |l^, the solid curve gives the branching fraction B{H — > 
Il^l) + B{H ^ again for fixed mu^ = 300 GeV and assuming three degenerate 

generations with masses m^-^ = 100 GeV and mj>^ determined by the relevant relations of 
Eq. (0), which is given in the next section. We see that the branching ratio is enhanced 
relative to the previous case, as the SU(2)l gauge couplings now contribute to the decay 
process, and decays to both charged sleptons and sneutrinos are now open. However, the 
branching ratios again drop rapidly for increasing tan/3. In Fig. |1^, we give results for 
tan (3 bounds with = 1 TeV, = 300 GeV, and all systematic errors included, and 
including the effects of SUSY decays to (a) right-handed sleptons with m^^ = 100 GeV, and 
(b) left-handed sleptons with mi = 100 GeV. As one can see, the loss in statistics is not 



very significant in both cases, and Fig. is almost unchanged. 

If decays to both left- and right-handed sleptons are allowed, the branching ratios of 
H^, A, and H are all altered. Decays to left-right pairs involve the /i parameter, as well as 
the trilinear scalar couplings. If these parameters are not measured, they may contribute 
large systematic errors to the measurement of tan (3. Of course, these parameters may also 
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be measured in different processes, for example, from chargino and neutralino masses for n 
and left-right mixings for the trilinear scalar coupling. A complete analysis would therefore 
require a simultaneous fit to all of these parameters. 

Finally, we briefly consider decays to charginos and neutralinos. These decays have been 
considered in detail , and have been shown to be dominant in some regions of parameter 
space. However, if only decays to the lighter two neutralinos and the lighter chargino are 
available, and these are either all gaugino-like or all Higgsino-like, as is often the case, these 
decays are suppressed by mixing angles. If we are in the mixed region, these decay rates 
may be large, but in this case, all six charginos and neutralinos should be produced, and 
the phenomenology is quite rich and complicated. 



VI. DETERMINING SOFT SCALAR MASSES 

As mentioned in Sec. |I|, the parameter tan f5 plays an important role in determining the 
masses and interactions of many supersymmetric particles. A measurement of tan f3 from the 
Higgs scalar sector is therefore valuable for constraining other supersymmetric parameters 
of the theory, or for testing SUSY relations in another sector. In this section, we present 
one simple example, namely, the determination of soft SUSY breaking scalar masses. 

The physical masses of sleptons and squarks are given (in first generation notation) by 

= m| + + m|(i - | sin^ 9w) cos 2/3 , 
= m| + ml + m|(— ^ + | sin^ 6w) cos 2/3 , 



L = m| + m^ + m|(|sin2 6'i^)cos2/5 

^2 _ ^2 I ^2 I ^2/ lo;^2, 



TTij = + + m|(— i sin^ ^^y) cos2/3 , (14) 
mj^ = mj^ + ml + m|(-| + sin^ 6w) cos 2/3 , 
ml^ = m| + m|(|) cos2/3 , 

= m|, + ml + m|(— sin^ 9w) cos 2/3 , 

where mg, m^, m^,, mi, and m^ are the soft SUSY breaking scalar masses. In these 
relations, possible mixings among sfermions are neglected. In fact, such mixings may be 
large and lead to a variety of new phenomena that may be probed at future e"'"e~ colliders. 
Left-right mixing, which may be large for third generation sfermions, has been analyzed for 
scalar taus [^, and intergenerational slepton mixing has also been studied recently pi] . 



For simplicity, however, we assume in this section that these effects are absent. 

As emphasized in Ref . P] , the pattern of soft SUSY breaking parameters is a reflection of 
the SUSY breaking mechanism, and so accurate determinations of the soft SUSY breaking 
masses may provide insights into the physics of SUSY breaking. In addition, accurate 
measurements of the sfermion masses may help determine the gauge and/or flavor structures 



at higher energies As can be seen in Eq. (p!4[) , an accurate measurement of the soft 

scalar masses requires precise measurements of both the physical sfermion masses and tan f3. 
If tan f3 is completely unknown, the uncertainty in the soft scalar mass is considerably greater 
than the expected uncertainty from measurements of the physical masses. For example, if 
mg^ = 100 GeV, the variation in for 1 < tan/? < 60 is 10 GeV. On the other hand, 
slepton and squark masses may be measured at e+e~ colliders without parameter unification 



assumptions with a fractional error of 1-2% |21,29,34]. One might hope, therefore, that a 
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measurement of tan [3 from the Higgs sector would reduce the uncertainty from tan /5 to a 
comparable level. 

In Fig. we plot contours of constant 

Am = |m^(tan/5) — m^(tan/3')| , (15) 

for fixed physical mass = 100 GeV. For other sfermion species and masses, 

the contour labels (0.5, 1, 2, and 3 GeV) should be multiplied by approximately 
(|F|/ sin^ 6'h/)(100 GeV/nij), where F is the appropriate function of sin^ in parenthe- 
ses in Eq. (|14D- (For example, for the right-handed slepton, F = — sin^6'iy.) The soft mass 
depends on tan/5 only through cos 2/5. The dependence is therefore very slight for large 
tan P, and, as may be seen in Fig. |2^, a bound such as tan /5 ^ 6 is already very powerful 
for the purposes of determining soft mass parameters. Comparing this with the bounds that 
may be achieved from Higgs scalars, we see that for tan/3 > 3, the uncertainty in m^. from 
tan/3 is reduced below that from the physical mass measurement. 

An independent measurement of tan j3 may also be possible in the scalar sector if the 
masses of both sfermions of a left-handed doublet are measured. (Such a measurement is 
not necessarily easy, even if sfermions are kinematically accessible — in the slepton sector, 
sneutrinos may decay invisibly; in the squark sector, such a measurement requires a de- 
termination of quark flavor.) In this case, a comparison of the two tan/3 determinations 



constitutes a highly model-independent test of SUSY |^ , without any assumptions of GUT 
or SUGRA relations. 



VII. IMPLICATIONS FOR TOP SQUARKS 

In this section, we discuss another application of the tan/3 measurement, namely, the 
application to the top squark sector. Top squarks are singled out by the large top Yukawa 
coupling, which implies that significant left-right stop mixing is generic, and that radiative 
corrections from top-stop loops are highly significant in determining the properties of the 
Higgs bosons. For these reasons, precise measurements of the parameters in the Higgs sector 
may allow us to constrain parameters in the stop sector. 

In the absence of left-right stop mixing, the leading radiative corrections to the CP-even 
Higgs sector were given in Eq. (^. In general, however, all parameters of the top squark 
sector enter. The top squark mass matrix is 

^2 / -|- -|- m|(| — I sin^ 6*1^) cos 2/3 mt{At — fi cot P) 

* Y mt(Aj — /icot/3) m|, m|(| sin^ ^vf) cos2/3 

where mg and mf are soft SUSY breaking parameters discussed in the previous section, fi 
is the supersymmetric Higgs mass parameter given in Eq. (^, and At is the top trilinear 
scalar coupling. The physical top squark masses are the eigenvalues of this matrix, and we 
denote the lighter and heavier top squarks as ti and t2, respectively. For low and moderate 
tan /3, where the bottom squark contributions may be neglected, the CP-even Higgs masses 
are then determined by ttia and tan/3 at tree level, and mg, ruf, fi, and At, all of which 
enter through radiative corrections. 
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In Fig. ^ we plot rrih as a function of m^^ for various values of tan (3 and At. Here, 
we fix /i = 200 GeV and = 300 GeV, and we take uiq = nif for simplicity. We see 
that there is strong dependence of m/j on the various top squark parameters, implying that 
we may be able to constrain new parameters with a measurement of tan/3. For example, 
if /i is measured from the gaugino-Higgsino sector, a measurement of tan/5, along with 
measurements of m^, rni^ and m/j, may allow us to constrain At, a parameter that otherwise 
may be rather difficult to measure without model-dependent assumptions. (Note that we 
have assumed the relation ttIq = rrif for simplicity. This relation may be relaxed if we 
also measure m^^ and impose this as an additional constraint.) On the other hand, from 
the figure, we also see an asymptotic behavior for large m^^ — if the soft SUSY breaking 
parameters mg and rrif dominate the left-right mixing terms, nih is simply a function of 
tan/?, niA, and m^^. Thus, assuming that this holds, even if top squarks are too heavy to 
be discovered at either a future e~^e~ collider or the LHC, we may be able to place upper 
bounds on their masses by measuring ruh, tua, and tan/5. It is important to note, however, 
that if At and fi may be arbitrarily large, one cannot draw such a conclusion. 

In this section, we have not considered quantitatively the results that may be achieved. 
Clearly, measurements of many parameters enter the analyses suggested here, and an overall 
fit to the relevant parameters will be necessary in a complete analysis. However, the exam- 
ple of the top squark sector illustrates at least qualitatively the possibility of applying an 
accurate tan (3 measurement to interesting determinations of parameters in other sectors of 
supersymmetric models. 



VIII. CONCLUSIONS 

In this study we have considered the prospects for measuring tan (3 through heavy Higgs 
scalar production and decay at a future e^e~ collider. The branching ratios of heavy Higgs 
scalars are strongly dependent on tan /?. In addition, we have seen that Higgs signals typically 
have many b quarks in the final state, which, given the excellent 6-tagging efficiencies and 
purities expected, allows them to be separated from SM backgrounds in a number of different 
channels. The cross sections from these channels allow us to significantly constrain the 
parameter space, and in particular, to place bounds on tan /3. 

We have relied on experimental measurements wherever possible. The neutral Higgs 
sector is subject to large radiative corrections, depending strongly, for example, on top squark 
masses and mixings. In our analysis, we have not assumed that such radiative corrections 
are small. Instead, we treat the Higgs scalar masses, the parameter cos^(/5 — a), and the 
vertex Xnhh as independent quantities, constrained only by experimental measurements. 
In addition, we have avoided assumptions of SUSY parameter unification. The analysis 
method is therefore formally applicable to models with arbitrary radiative corrections. We 
have, however, assumed a minimal Higgs sector throughout this analysis. If additional Higgs 
fields are present, determinations of tan/5 from the various channels may not be consistent, 
and the effects of a non-minimal Higgs sector may also be detected by an analysis similar 



to this one 35 



We have considered scenarios with a variety of beam energies, Higgs masses, integrated 
luminosities, 6-tagging efficiencies, and assumptions about systematic errors, and have also 
considered the impact of SUSY decay modes. For all scenarios considered, we find that the 
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strong dependence of heavy Higgs branching ratios on tan (3 allows stringent constraints for 
moderate tan/3. These results imply that for tan/? ^ 3, the soft scalar mass parameter 
determinations are likely to be limited by the precision of the corresponding physical scalar 
mass measurements. In addition, we have seen that the three body cross section tbH^ 
grows rapidly for large tan j3, and is large enough for some scenarios to provide interesting 
constraints for large tan (3. These bounds allow one to confirm or exclude Yukawa unification 
assumptions. 

The parameter tan jS may also be constrained by other processes. For example, for low 
tan/5, chargino production at a linear e~^e~ collider may provide stringent constraints on 
tan/5 This requires a sufficient Higgsino component in the chargino, and the bound 

deteriorates for moderate and high values of tan/?. For high tan/?, there are a number of 
possible probes. If staus are pair produced at a future linear collider, tl-tr mixing may 
be able to measure tan j3 sensitively in a range determined by the Higgsino component of the 
lightest neutralino. Alternatively, the magnetic dipole moment of the muon, {g — 2)^, may 
be sensitive to tan/5 ^ 20 for slepton masses m^- < 300 GeV [^. In addition, the discovery 



oi H,A ^ TT at the LHC may be used to set the lower bound tan/5 ^ 10 [^. In general, 
however, the heavy Higgs sector appears to be one of the most challenging for the LHC, and 
the discovery and study of heavy Higgs bosons there may be difficult Of course, in the 
ideal case that measurements confirm a particularly simple model, for example, the so-called 
minimal supergravity model, in which all supersymmetric particle masses and interactions 
are determined by only five additional parameters, studies have shown that highly accurate 



measurements of tan (3 may be obtained both at the LHC and NLC . 

It is interesting to note, however, that the power of these other methods is usually 
greatest for low or high values of tan /5. In addition, these measurements all involve many 
other SUSY parameters and require certain conditions to be applicable. In contrast, the 
heavy Higgs measurement is most sensitive in the range where these other measurements 
are weak, and is relatively free of other assumptions. It is clear, however, that no one process 
is powerful throughout the range of tan (3 and for all scenarios. Of course, if more than one 
test is available, their consistency will be an important test of SUSY. 
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APPENDIX: DECAY WIDTHS 



In this appendix we give formulas for the decay widths of the heavy Higgs scalars 
to quarks and leptons for reference. Additional decay modes are H~ W~h, H — »• 
W~^W~ , ZZ, hh, AA, and A Zh, as well as SUSY decay modes involving squarks, slep- 
tons, charginos and neutralinos. All of these decay widths may be found in Appendix B of 



Ref. Ill 



The charged Higgs boson decay width for fermion pairs fufd is 



m 



fu 



m 



(Al) 



where Nc is the number of color, and we have approximated ruf^ <^ in the phase space 
factor. 

For A and H decays to //, the width is given by 
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where the coefficient C and exponent p are specified as follows: 
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TABLES 



TABLE L The cut parameters for our simulation. 

^/s = 500 GeV ^ = 1 TeV 



AE^± 30 GeV 30 GeV 

AE+± 5 GeV 10 GeV 

AE^^ 15 GeV 25 GeV 

A£;2b 20 GeV 35 GeV 

Arriz 30 GeV 40 GeV 

Am| 5 GeV 10 GeV 

Amw 10 GeV 15 GeV 



TABLE II. The standard model background for channels 1 - 8 (in fb). Here, we assume 
eb = 60% and Cc = 2.6%. 





1 


2 


3 


4 


5 


6 


7 


8 


^/s = 500 GeV 


0.11 


0.071 


1.5 


6.2 


1.2 


0.12 


0.17 


0.0053 


= 1 TeV 


0.36 


0.35 


1.0 


3.0 


0.47 


0.22 


0.30 


0.0095 



TABLE III. Cumulative efficiency of cuts la - le for = 500 GeV. Here, we take m^± = 200 
GeV and tan/3 = 5. The 6-tagging efficiency is not included. 



cut H+H- tt AH 

Ta 090 1.2 X 10-3 4.5 x lO'^ 

lb 0.82 1.1 X 10-3 4.5 X 10-^ 

Ic 0.51 5.7 X 10"^ 3.8 x 10"^ 

Id 0.51 5.7 X 10-4 3.8 x 10-^ 

le 0.45 5.6 x 10-^ 3.4 x 10-^ 



TABLE IV. Cumulative efficiency of cuts 2a - 2e for -s/s = 500 GeV. Here, we take mjj± = 200 
GeV and tan (3 = 60. The tbH"^ signal is normalized to 1 after cut 2a. The 6-tagging efficiency is 
not included. 



cut tbH^ tt AH 

'2a LOO 095 033 

2b 0.65 3.1 X 10-4 0.27 

2c 0.38 2.2 X 10-4 0.27 

2d 0.38 2.1 X 10-4 0.26 

2e 0.33 1.5 x 10-4 0.24 
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FIGURES 




500 1000 1500 



^[GeV] 



FIG. 1. Production cross sections for e^e~ H^H^ (solid) and AH (dashed) for mjj± = 200, 
300, and 400 GeV from left to right, and Zh (dotted) for = 200 GeV. (The Zh cross sections 
for mjj± = 300 and 400 GeV are virtually identical.) The 1-loop radiative correction given in 
Eqs. (^ and (^) has been included with = 1 TeV, and we have set tan /3 = 5; the dependence 
on tan (3 is very weak for mu± ^ 200 GeV. 
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FIG. 2. The three Feynman diagrams contributing to the three body final state tbH . Three 
similar diagrams contribute to tbH~^. 
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FIG. 3. The thrcc-body cross section cr(e+e tbH ) + o{e'^e thH^) for m//± = 200 
GeV (solid), 300 GeV (dotted), and 400 GeV (dashed), with tan/3 = 60. We require 
Et + Eb > 1.02y^/2 to separate this mode from the two-body production of H^H~ followed 
by tb. 
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FIG. 4. The branching ratios for H'^ to rP (sohd), tb (dashed), and W~h (dotted) for (a) 
mH± = 200 GeV and (b) mfj± = 400 GeV. The leading radiative correction of Eqs. (^) and (Q) 
is included with = 1 TeV in calculating the remaining Higgs masses and mixings. The running 
quark mass = 3.2 GeV has been used, and all SUSY decay modes are assumed suppressed. 
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FIG. 5. Same as in Fig. ^, but for CP-odd scalar A decays to rr (solid), bb (short dash), tt 
(long dash), and Zh (dotted). 
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FIG. 6. Same as in Fig. but for neutral Higgs scalar H decays to rr (solid), bb (short dash), 
tt (long dash), hh (dotted), Ty"'"W^^ (dot - short dash), and ZZ (dot - long dash). 
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FIG. 7. Schematic pictures of signal events H'^H tbrv and tbH thru, and background 
tt bqq'bW. The crosses mark the interaction point. 
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FIG. 8. Signal cross sections for mH± = 200 GeV, tan/3 = 5, and -s/s = 500 GeV, in each 

channel: (1) 2b + I + q's + cuts la-le ("ff+if"" channel) (solid), (2) 2b + I + q's + cuts 2a-2e 
{HbH^" channel) (dotted), (3) 3b +11 (+g's) (short dash), (4) 3b + 0,2, 3,... I (+g's) (long dash), 
(5) 46 (dot-short dash), 6) 46 + II (+g's) (dot-long dash), (7) 46 + 0, 2, 3, . . . I (+q's) (but not 46) 
(short dash-long dash), and (8) 5b {+ I + g's) (dash-dot-dot). All branching fractions, tagging 
efficiencies, and cuts have been included. 
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FIG. 9. Signal cross sections as in Fig. ^, but for m/^± = 200 GeV and ^/s = 1 TeV. The cross 
section for channel 8 is less than 10~^ fb, and therefore does not appear. 
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FIG. 11. Signal cross sections as in Fig. but for mff± = 400 GeV and ^/s = 1 TeV. 
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FIG. 12. 95% C.L. bounds on tan/3 for y/^ = 500 GeV, mjj± = 200 GeV, Sb = 60%, and four 
integrated luminosities: 25, 50, 100, and 200 fb~^ (from outside to inside). For a fixed underlying 
value of tan /?, the values of tan (3' determined by the appropriate contours are the upper and lower 
bounds that may be set experimentally. 
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FIG. 17. 95% C.L. bounds on tan /3 for (a) mj^± = 300 GeV and (b) mH± = 400 GeV, ^/s = 1 
TeV, ej, = 60%, and four integrated luminosities: 100, 200, 400 and 800 fb~^ with all systematic 
uncertainties omitted. 
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FIG. 18. Branching ratios B{H 1*Jr) (solid) and B{H llli) + B{H vIDl) (dashed). 
The H mass and couphngs are determined for fixed mjj± = 300 GeV, and include the radiative 
correction of Eqs. ^ and @) with rrij = 1 TeV. For the solid (dashed) curve, the three generations 
of right- (left-) handed sleptons are assumed degenerate with mass mj = 100 GeV ("^-f^ = 100 
GeV), and all other sparticles decay modes are assumed closed. 
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FIG. 19. 95% C.L. bounds on tan/? for ^ = 1 TeV, mj^± = 300 GeV, eb = 60%, and four 
integrated luminosities: 100, 200, 400, and 800 fb~^, with (a) only the right-handed slepton decays 
open with m^~^ = 100, and (b) only the left-handed slepton decays open with m^-^ = 100 GeV. All 
systematic errors are included. 
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FIG. 20. Contours of constant Am = 0.5 GeV (solid), 1 GeV (dotted), 2 GeV (short-dashed), 
and 3 GeV (long-dashed), where Am = |mg(tan/3) — m^(tan/3')|, for fixed physical mass 
me„ = 100 GeV. 
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FIG. 21. The lightest Higgs mass m/j as a function of the mass of the Hghter stop m^^ for 
tan/3 = 3, 5, and 10. Here, we fix = 200 GeV, itia = 300 GeV, and At = (sohd hnes), 
At = 200 GeV (dashed hnes), and At = -200 GeV (dotted hnes). 
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